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APPLICATION OF TEE MATCHED CONIC MODEL 
IN THE STUDY OF 
CIRCUMLUNAR TRAJECTORIES 


ABSTRACT 


Matched conic and integrated circumlunar trajectories with the 
tame constraints are compared at injection, pericynthion and perigee. 
Results indicate that the matched conic model may he used for design 
purposes and for obtaining trends of circumlunar trajectories. In 
addition, an iteration scheme is presented to obtain any desired 
integrated trajectory. The scheme utilizes the matched conic model 
in conjunction with the integrated n-body model. Convergence is rapid. 


INTRODUCTION 


Two problems encountered in the study of circumlunar trajectories 
are: (l) obtaining a reliable analytic model of the earth-moon system 

and (2) obtaining a fast converging iteration scheme for the determi- 
nation of initial conditions for precision trajectories. These problems 
arise because there is no known analytical solution to the general 
n-body problem where there are more than two bodies involved. It is 
therefore desirable to obtain an approximate earth-moon model which 
brs an analytical solution and yields results which are good 
approximations to precision (integrated) results. This model 
may then be used for parametric studies, as an initial condition 
generator, and to form a basis for an iteration scheme to converge 
upon initial conditions for a more accurate model. The purpose of 
this paper is to show that the matched conic is a model which closely 
approximates the results of precision trajectories and to show how 
the matched conic model may be linked with a precision model to give 
a fast converging, reliable iteration scheme. 

These applications divide the discussion into two sections, 
both of which refer to the matched conic and 7-tody trajectory 
models. A brief discussion of each is given to familiarize the 
reader with the models. For more detailed information concerning 
the models, see references 1 and 2. 
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SYMBOLS 


N 

R 

S 

T 

V 
h 
i 
r 

n 

Y 
0 

Y 

U) 

Prefixes 

A 

V 

Subscripts 

c 

d 


trans lunar or transearth trajectory above 
the moon orbital plane, deg 

distance of closest approach (r in example) 

translunar or transearth trajectory below 
the moon orbital plene, deg 

time of flight, hr 

velocity, ft/sec 

injection altitude, ft 

inclination 

radius 

angle of ascending node, deg 
flight, path angle, deg 
lead angle, deg 
argument of radius, deg 
argument of periapsis, deg 

correction term for matched conic initial 
conditions 

error in integrated end conditions 

matched conic 
desired 


e 


earth 
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i integrated 

m (e. c. ) moon 

r return 

tl trans lunar 

Superscripts : 

* constrained element 

DISCUSSION 


Matched Conic 

The particular matched conic model used in this discussion 
was initially developed hy the Martin Company for MSA in con- 
nection with the Apollo study contract (ref. l). The model 
simplifies the complicated earth-moon system to just two point 
masses which represent the earth and moon. The motion of a 
third body (particle) in this system is governed completely by 
the gravitational force of either the earth or moon. To determine 
which force is to be used, a fictitious sphere, called the moon's 
sphere of influence, is constructed such that it encloses the 
moon and moves with it. The radius of the sphere is constant and 
was chosen empirically by comparison with integrated trajectories. 

If the third body is inside the sphere of influence, the moon's 
gravitational force is used; otherwise, the earth's gravity is used. 
Since the force acting on the third body is always inverse square, 
the orbits are always solutions to the Kepler problem or conic 
sections. Thus, a lunar trajectory in this system is a series of 
conics which are matched at the sphere of influence. Specifically 
for circumlunar trajectories, the particle would travel out to the 
sphere of influence on an earth- focused ellipse determined by the 
initial conditions. At the sphere, the velocity and position vectors 
are transformed to the moving moon system which results in a hyper- 
bolic orbit about the moon. When the particle again reaches the 
sphere of influence its orbit is transformed back to the earth 
reference. (See fig. 1.) 

All the orbital parameters are confuted by the usual conic 
formulas with the exception of the velocity at pericynthion 
which is computed by Jacobi's Integral to the Restricted Three- 
Body Problem in order to give a better approximation to integrated 
results . 
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Although analytic expression exists for a particle's trajectory 
in this model, the expressions relating initial and final conditions 
are transcendental and require iteration for solution. However, 
because of the simplicity of the equations, the iteration can be 
performed very quickly on a digital computer. A solution where 
three end conditions are satisfied requires about 10 seconds on the 
IBM 7090. 


7-Body 

The 7-hody model used for comparison is the NASA Inter- 
planetary program (ref. 2). Briefly this model computes a 
particle trajectory by Encke's Method subject to the gravita- 
tional forces of the following bodies: Earth, Moon, Sum, Venus, 

Mars, and Jupiter. Perturbations due to the non-spherical 
gravitational fields of the earth and moon are also taken into 
account. Initial conditions are given to start the program and 
a trajectory results upon numerical integration. The computa- 
tional time for a single trajectory of 150 hours requires about 
1.5 minutes on the IBM 7090. 

Results of Trajectory Comparisons 

A comparison was made with circumlunar trajectories which were 
computed on the matched conic and 7-hody models that satisfied the 
same constraints. The trajectories are the figure eight type as 
shown in figure 2. The constraints were imposed at three points 
of particular interest - at the injection, at pericynthion, and 
at the return perigee. At injection, there are three constant 
constraints imposed: radius (r), flight-path angle (y), and 

translunar inclination (i^.) . At pericynthion, one constraint, 

a constraint, a constant pericynthion radius (r ffl ), is imposed. 

At perigee, two constant constraints, perigee (r ) and return 

inclination (i re ), are imposed. The geometry and definition of 

the constraints at injection, pericynthion, and perigee are given 
in figures 5 to 5, respectively. Specifying the six constraints 
and the initial time determines a circumlunar trajectory uniquely 
on both models. A total of 18 elements are defined at the three 
points, six at each point. Six of the 18 elements have the same 
constant values on both models because of the constraints; however, 
the remaining 12 elements will generally have different values on 
the two models. The 12 unconstrained elements are called free 
elements. At injection, the free elements are velocity (V), lead 
angle (0), and argument of injection (Y). At pericynthion, the 
five free elements are time from injection (T ), velocity (V ), 
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inclination (i^), angle of ascending node (n^), and argument of 

pericynthion The four free elements at perigee are time 

from injection (T m ), velocity (V g ), angle of ascending node 

and argument of perigee (u) e ). The deviations of the matched conic 

elements from the 7-body elements will give a good indication of 
the accuracy of the matched conic model. The definition of the 
free elements at the three points is defined and denoted on 
figures 3 to 5- A summary of the free and constrained elements 
is given in table I. 

A sample matched conic and 7-body comparison is given in 
table II. The injection time for all the trajectories presented 
is Greenwich midnight January 15, 1967. The moon is near apogee 
at the time of arrival at pericynthion. The two trajectories have 
a pericynthion altitude of 100 nautical miles and a perigee altitude 
of 20 nautical miles. The inclination at perigee is 5° south. The 
differences are computed by subtracting the matched conic element 
from the corresponding 7-body element. At injection, the differences 
in the constrained elements are 0 because they are direct inputs. 

Small difference in the constraints exists at pericynthion and perigee 
because of the small convergence tolerances in the iteration scheme. 

To give further evidence that the differences are small, 

trajectories were run on the two models where one of the six 

constraints was varied. The return inclination (i ) was varied 

' re 

between 6o°S and 60°N. The trajectories were injected at an 
altitude of 600,000 ft, 0° flight-path angle (y), and 1°N 
translunar inclination (i ). The pericynthion altitude was 

100 nautical miles and the perigee altitude 20 nautical miles. 

Plots of the free elements as functions of i for both models 

re 

at injection, pericynthion and perigee are given in figures 6 to 8, 
respectively. It is seen from the plots that the shapes of the 
matched conic curves are very near the same shapes as the corre- 
sponding precision curves. In addition, the matched conic curves 
are displaced from the integrated curves by small increments which 
are nearly constant. Hence a small constant correction term may 
be added to the matched conic value to give even a better approxi- 
mation. These facts are the foundation of the iteration scheme 
presented in the following section. 

In addition to the few comparisons presented here, the author 
has made many spot checks with variations in the six constraints. 

In all cases, the matched conic performed as well as presented here. 
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On this basis, it seems justified to infer that the differences 
remain small for all figure eight circumlunar trajectories. 

Iteration Scheme 

The problem of calculating a trajectory which satisfies certain 
constraints along its path is called a boundary- value problem. The 
solution to the circumlunar boundary- value problem requires a 
numerical solution because of the absence of an analytic solution 
to the n-body problem. Generally, in obtaining a numerical solution, 
a first guess is made to the initial conditions and a trajectory is 
computed by numerical integration. The values of the end conditions 
are compared with the desired ones. If any of the end conditions 
are not within the allowable error, another guess of the initial 
conditions is made. The procedure of guessing (iteration) and 
integrating continues until the desired end conditions are met or 
until the results indicate that the desired solution does not exist. 
Obviously the rate of convergence to the proper initial conditions 
depends upon the particular iteration scheme. Since the computa- 
tional time of a single 7-body integrated trajectory is about 

1^ minutes, an iteration scheme is desired so that only a few 

integrated trajectories need be computed for convergence. 

Ordinarily, an iterative scheme would rely on a near-linear 
relationship existing between the initial conditions and end 
conditions. However, for eircimlunar trajectories it is well known 
that the relationship Is very nonlinear. 

The A iteration scheme presented here relies on the near 
constant differences existing between matched conic and 7-body 
trajectories which were presented in the previous section of this 
paper. This means that a set of constant correction terms may be 
used to update the matched conic initial conditions. Also, changes 
in the initial conditions will affect the end condition in the same 
manner for both models. 

The simple theory of the A iteration scheme can be^t be shown 
by figure 9 which exhibits a single variable iteration. The figure 
shows two curves of injection velocity as a function of return 
inclination. The correction term AV is assumed to be known from 
a previous trajectory comparison and is used throughout the iter- 
ation. For the first guess, the matched conic velocity is computed. 
This matched conic velocity corresponds to point 1 q on the matched 

conic curve. The term AV is added to the matched conic velocity, 
and this new velocity is used in the integrating program. Upon 
integration, the trajectory arrives at the point 1^ on the 
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integrated curve. An error now exists in the return inclination. 

Since trends are nearly the same, e change of V 1 in the return 

inclination of the matched conic will produce nearly the same 
change in the return inclination of the integrated. Therefore, 
the inclination is moved to the right by on the matched conic 

curve, which gives point 2 . Again TV is added and this gives 

point 2^ on the integrated curve. A much smaller error Vg now 

exists. The matched conic is changed by Vg and the procedure 

continues as before until V is within the allowable tolerance. 

(Note that the error in the integrated end condition is always 
added algebraically to the matched conic end condition. ) If the 
AV used was the exact one for the desired inclination, no iteration 
would be necessary. But, in general, the AV used will not be the 
exact one. It Is immediately evident that the function of AV is 
to make the iteration occur between closer curves. For multiple 
convergence, the procedure is the same; correction terms are added 
to the free initial conditions, and the errors are added to the 
end conditions of the matched conic. In general, the correction 
terms are necessary because uncorrected initial conditions give 
rise to integrated trajectories which are usually grossly in error. 

A step-by-step procedure for a multiple convergence follows . 

For clarity, the method will obtain initial conditions for a 
circumlunar trajectory which has desired value of pericynthion, 
perigee, and return inclination (see fig. 2). Injection is from 
a fixed translunar inclination, altitude, and flight-path angle. 

The injection geometry is shown in figure 3. A set of correction 
A's is known from a previous trajectory comparison and is used 
throughout the iteration. In case a set of A's is not available 
or a closer set is desired, they may be obtained by the following 
steps: 

Step 1 : Run a matched conic trajectory with the desired 

constraints . 

Step 2 : Run integrated trajectory using matched conic 

initial conditions. 

Step 3: Vary one of the integrated free initial conditions 

(V, ¥, or 0) until a trajectory has the proper 
return perigee (no more than five integrated 
trajectories are required in using linear scheme). 
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Step 4 : Run a. matched conic vhich has the same pericynthion 

and return inclination as the integrated trajectory 
with the proper perigee. 

Step 5 ? Compute the difference in the free initial conditions: 

AV 
AY 

A0 

Die same set of A's has been adequate for all the free return 
trajectories computed by the author to date. However, for faster 
convergence, the A's shown in figure 10 should be used. With a 
set of A's available, the steps in the A iteration scheme are 
outlined below and summarized in block diagram form in figure 11. 

S tep 1 : Run matched conic with desired end conditions. Die 

output is a set of initial conditions. 

Step 2 ; Correct initial conditions by A’s. 

Step 3 : Compute integrated trajectory by using corrected 

initial conditions. 

Step 4 : Compute the error differences (A's) in the end 

conditions (desired value minus integrated value). 

Step 5 : If errors are small enough, iteration is terminated. 

If not, errors are added to previous matched conic 
end conditions. 

Step 6 : Repeat the process, but in Step 1 use new matched 

conic end conditions. 

In general, the number of iterations necessary depends upon the al- 
lowable end condition tolerances and accuracy of the A's. Ordinari- 
ly 6 or 7 iterations are required (10 minutes computer time). 

To illustrate, a sample iteration is shown in figure 12. Die 
constraints are at the top of the figure together with the cor- 
rection A's. The diagram is largely self-explanatory, however, 
the following explanation is given for the convergence to return 
perigee radius. The other parameters, pericynthion radius and 
return inclination, are iterated simultaneously with the perigee 
radius and in the same fashion. For the first iteration, the 




9 


desired value of 3, 464 nautical miles was used in the m tched conic . 
The resulting integrated value was 3*164 nautical miles or 300 
nautical miles too low. Hence, for the next iteration, perigee on 
the matched conic was increased by 300 nautical miles to equal 
3,764 nautical miles. Instead of increasing by 300 nautical miles, 
the integrated value only increased 274 nautical miles or, still an 
error of 26 nautical miles. On the next iteration the matched conic 
was further increased by 26 nautical miles to a perigee of 
3,790 nautical miles. The use of a corrected matched conic initial 
condition for this trajectory resulted in an integrated trajectory 
with an error in perigee of only 1.1 nautical miles. In each 
successive iteration the errors were reduced. Die magnitude of 
the errors shows that perigee radius is the most sensitive. The 
exact A's for the trajectory are obtained by subtracting the first 
matched conic from the last corrected matched conic. Hiese A's 
are shown at the bottom of the figure. 


CCMCLDDIHG REMAKES 


Hie matched conic elements approximate precision elements 
adequately for many types of parametric studies. An even better 
approximation may be obtained by adding constant correction terms 
to the conic element. Hie A iteration scheme has been tested and 
in all cases proved to be superior to linear perturbation type 
iteration schemes in both dependability and computational time. 


REFERENCES 


1. Anon: Apollo Final Report - Trajectory Analysis. Martin Rep. 

No. ER 12003, June 196 I. 

Fines, Samuel, an- 3 Wolf, Henry: Interplanetary Trajectory 

by the Ehcke Method Programmed for the IBM 704 and 7090. ~* 

Republic Aviation Corporation Rep. No. RAC- 656 - 451 , 

Dec. i 960 . 


2 . 



TABLT I.- GflIITAL ELEMENTS 


Point 


Elements 


Fixed constrained 


Free 


Injection 

Pericynthton 


r, r, i 

r m 

V ‘re 


tl 


v, Hr, 0 


T , i , V , n , o> 
m' m m' m' m 


V V n e- 


Return perigee 









Injection 


TABLE II.- TRAJECTORY COMPARISON 


Point 




Element 

Matched conic 

7-Body 

Difference 

V, ft/sec 

35,956.901 

35,959-242 

2.341 

V, deg 

7-75578 

7.330060 

-.42572 

<f, deg 

38.14462 

38.97043 

.82581 

*tl*' *** 

30 N 

30 N 

0 

y*, deg 

0 

0 

0 

r*, ft 

21,525,244 

21,525,244 

0 























TABLE II.- TRAJECTORY COMPARISON - Concluded 


Point 

Element 

Matched conic 

7-Body 

Difference 


V' ft 

21, 133,527 

21,133,948 

421 


i re *, deg 

5.186 S. 

5.184 s. 

.003 

a> 

<t> 

T e , hr 

148.85 

147.98 

-.85 

t 

£ 

O e , deg 

177-38 

176.87 

-.51 


a> e , deg 

-II.78 

-10.98 

-.79 


V , ft/sec 

56, 626 

36,631 

5 
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Figure 2.~~ Free return circumlunar trajectory 




geometry 
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Figure 4. “ Pericynthion parameters 




Note: XYZ axes same as XYZ for injection 


Figure 5. - Return parameters at perigee 
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Figure 7 - Continued 
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Figure 7 . - Continued 








Injection velocity 



Desired inclination 


Return inclination 


Figure 9. - Single variable A iteration 


AV, ft/sec 






Figure 11.- Flow diagram of A interation schema 




Injection time, Jan. 15, 1967 AV = 2.30 ft/sec 

h = Injection altitude = 600,000 ft - 0.38° 

Translunar inclination 20° N q 730 

Y-- 0 

Desired pericynthion (r m ) - 1,038.5 nautical miles (100 nautical miles altitude) 

Desired perigee (r e ) = 3,464.0 nautical miles (20 nautical miles altitude) 

Desired return inclination = 5° S 


First iteration Second iteration Third iteration 




Errors 

Errors 


Errors 

V 

r m -16 mi 

V r m 2mi 

V 

r m ”0.4 

V 

r e 300 mi 

V 26mi 

V 

r e l.lmi 

V 

i r 1 . 2 ° 

V i r 0.6° 

V 

i r -.06° 


Actual A’s: AV = 2.257 ttM.3925 A# = 0.7859 


Figure 12.- Sample iteration 
















































